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INTRODUCTION 


Scisiiiii  wave  al  ti-nuat ion  in  rocks  is  a  very  important  i)ara!n(.'tcr  not  only  because  it  atfrets 
the  seismic  wave  propagation  through  the  earth,  but  also  la'cause  it  is  a  sensitive  indicator 
of  rock  properties  under  various  conditions.  The  main  juirpose  of  measuring  attenuation  in 
tlie  laboratory  is  to  use  the  data  to  infer  rock  properties  at  vi-situ  scales  aiid  seismic  fi ■•quen- 
cics.  Unfortunately,  accurate  measurements  of  intrinsi<-  attenuation  are  difficult  to  oiTj.jn 
ii.  both  the  laboratory  and  fu-ld,  because  other  effects,  like  geometrical  spreading,  i.'/.-.r.dary 
t(  flections,  also  affect  waveform  amplitudes.  In  the  laboratory,  attenuation  is  usually  m.ca 
>ured  using  ultrasonic  pulse  propagation  and  resonant  bar  techniques.  For  the  propagation 
of  stress  pulses  ibrough  rock  samples,  attenuation  can  be  estimated  from  the  rise-time  of  the 
pulse  on-set  (Blair,  1982)  or  from  using  the  spectral  ratio  technique  (Toksbz  et  al..  1979). 
However,  as  pointed  out  by  Liu  (1988),  the  rise  time  is  strongly  affected  by  the  source  time 
function.  The  conversion  from  rise  time  to  attenuation  also  assumes  a  specific  attenuation 
model  (e.g.,  constant-Q  model,  Kjartansson,  1979),  which  may  become  invalid  in  the  pres¬ 
ence  of  fluid  saturation  (Jones,  1986).  The  spectral  ratio  method  derives  attenuation  from 
the  slope  of  a  line  fitted  to  the  logarithmic  ratio  of  the  spectra  of  two  waveforms.  Because 
waveforms  may  be  affected  by  reflections  from  sample  boundaries,  high  frequency  pulses  (~ 
1  MHz)  and  samples  with  large  lateral  dimensions  arc  used  in  the  measurements  (Toksoz  et. 
al.,  1979).  Even  so,  signals  may  still  have  to  be  truncated  to  remove  extra  arrivals  which 
may  significantly  affect  the  wave  spectra.  Moreover,  the  diffraction  effects  of  the  transducer 
source  may  nerxl  to  be  corrected  for  in  such  measurements  (  Lang  et  al.,  1990).  The  resonant 
l>ar  technique  utilizes  the  resonance  of  a  vibrating  cylindrical  rod.  The  resonant  frequency 
of  the  fundamental  mode  usually  occurs  between  3~10  kHz.  depending  on  the  length  of  the 
rod.  Because  of  the  very  different  frequency  ranges  of  the  ultrasonic  and  resonant  techniques, 
attenuations  obtained  from  the  two  techniques  can  be  significantly  different  (Blair.  1990). 
In  addition,  the  resonant  bar  me<isureinent  is  difficult  to  perform  under  pressure,  whereas 


I  hr  [.ul-  c  iiH‘lh.-<l  is  most,  suili-.l  loi  iiv»-  in  |,i.  smu.'  ||  js  <|.  s||,,li|r  to 

l:a\<:  a  Uchni(iue  tliat  iiK.-a.siiics  att ciiuaf ion  in  tin-  low  to  iiR-iiiniu  lic<in<  iuy  langr  in  ooi.  i 
to  nndcistand  attenuation  meclianisins  in  rocks  as  a  function  of  wavelengths  and  of  ia-s./a 
<  onditions.  Tlie  purpose  of  this  study  is  to  develop  a  technique  that  combines  the  advantage-.^ 
of  the  two  methods  and  opcratc*s  in  the  low  to  medium  frequency  range's. 

.■\  .lircct  way  of  reducing  the  measurement  frequency  in  the  pulse  proi)agation  mel!u/d  is 
b\  udiig  low  frequency  sources  and  increasing  the  propagation  length  between  source  and 
icci  i'.cr  transducers.  However,  bc-eausc  of  the  beam  spreading  of  the  transducer  radic.tic.n. 
ihc  lateral  dimension  of  the  sainjrle  must  be  increasc-rl  proportionally  to  avoid  the  contam¬ 
ination  of  (lie  direct  signal  by  the  reflections  from  lateral  boundaries.  This  re.suhs  in  a 
large  sample  volume  that  is  impractical  for  laboratory  measurements.  'I'his  problem  can  bi' 
a'.oi'led  bv  using  the  c.ylinder-shape<l  sample  as  used  in  the  re.sonaut  bar  metliod.  'I’he  u-( 
of  the  bar  geometry  in  the  pulse  propagation  method  has  two  major  advantages.  First,  the 
sanjple  length  can  be  chosen  according  the  wavelength  without  having  to  increase  lateral 
dimensions.  Thus  it  is  suited  for  use  in  pressure  vessels  and  with  saturated  samples.  Second, 
the  fundamental  mode  in  a  cylindrical  bar  is  a  low  frequency  wave  phenomenon  whose  prop¬ 
agation  and  dispersion  characteristics  are  well  understood  and  can  be  accurately  modeled. 
Because  the  dispersion  effects  of  the  waveguide  can  distort  the  source  signal  into  a  lung  wave 
train  (particularly  at  high  frequencies),  the  spectral  ratio  technique  is  not  suitable  for  thii. 
application.  This  method  is  most  accurate  with  short  duration  signals.  To  overcome  tliis. 
we  have  developed  a  waveform  inversion  technique  that  needs  only  the  first  few  cycles  uf  the 
waveforms  to  obtain  reliable  estimates  of  attenuation. 

In  the  following  studies,  the  propagation  characteristics  of  the  fundamental  wave  mi-dc  in 
.1  ■:}  lindrical  l^ar  will  be  discussed.  Then  the  procedure  for  the  wa\eform  imersion  technique 
wii!  be  formulated.  Finally,  the  procedure  is  applied  to  measme  attenuation  in  a  I’N'C 
materia]  and  in  Sierra  While  granite.  The  results  for  the  granite  arc  compared  with  those 
obtained  using  other  ttrchniquc's. 
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PROPAGATION  OF  THE  FUNDAMENTAL  WAVE  MODE  IN 

A  CYLINDRICAL  BAR 


III  a  lliin  cyliii<irical  rod  consisting  of  an  clastic  material,  jn-opagation  of  cxtcnsio.'-.a'i  •Aa.’.es 
goserned  by  the  Poclihaminci  ixpiation  (Kolskv.  Ib'i3) 

['2i/<i){in~  d-  k")Ji{la)Ji(ma)  —  {ni^  —  Jo{la)Ji{ina)  -  ik"hnJi{la)Jc{ma)  =  0  .  l  l) 

\\lii’rc  (I  is  bar  radius,  J„  (n  =  0,  1)  is  the  nth  order  Hesscl  function,  k  ;s  the  e.\tc:is;o;  al 
\va\'enumljcr,  and 

I  =  and  m  = 

are  radial  coinpressional  and  shear  wavenumbers,  and  are  compressional  and  shear 
'.elocities  re.spectively,  and  u  is  angular  frequencj'.  Although  Eq.  (1)  gives  rise  to  a  number 
of  (wtensional  wave  modes  in  the  bar  (Kolsky,  1953),  the  fundamental  mode  is  of  the  m.ost 
intiTcst  to  us.  At  very  low  frequencies  (a;  —►  0),  the  velocity  of  this  mode  approaches 
I  ^  =  \J E/ p,  where  E  is  the  Young’s  modulus  and  p  is  the  density  of  the  bar.  With  increasing 
frequency,  the  phase  velocity  of  the  mode  decreases.  Because  of  the  change  of  velocity  with 
frc<iuency  (i.e.,  velocity  dispersion),  the  waveform  of  the  mode  will  be  distorted  as  the  wave 
propagates  along  the  bar.  To  demonstrate  this  effect,  we  consider  the  spectrum  of  the 
propagating  wave  at  the  distance  x  away  from  the  source 

lE(w,  i)  =  5(w)  exp(zA-.r)  ,  l'2) 

wlii'if  .‘'’(u,’)  is  the  spectrum  of  the  transducer  source.  After  solving  Eq.  (1).  Eq.  (2)  (an 
lie  transformed  into  time  domain  to  obtain  waveforms  at  \arious  distances  i.  Figure  la 
■allows  the  phase  velocity  (dashed  curve  marked  'pha,«e'.  obtained  as  w/^')  gioup  \elocity 
(dashed  curve  marked  ‘group’,  obtained  as  du^/dk),  as  well  as  the  amplitude  of  tlie  souice 
spi  (  trum  S(u:)  (solid  cur\e)  as  functions  of  frequency  for  this  example.  Fhe  results  arc  calcu- 
l.il'  <1  for  a  bar  of  1  cm  radius  with  V),  •  Vll'l  m/s  anri  1,  --  J.aOO  '1  he  souice  is  a  Ri(  l'er 
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sourer-  (Ricker,  1953)  centered  around  70  kHz.  l-'igure  11)  sliows  (lie  .syntlictic  waveform^  ai 
the  distances  a.=0.  10,  20,  and  30  cm  from  the  source,  liecause  tlie  wave  spectrum  covers 
the  frequenc\'  range  where  significant  vcloritj'  dispersion  occurs,  tlic  waveform  is  gradu.iliy 
distorted  into  a  long  wave  train  as  it  propagates?  along  the  bar.  In  addition,  as  indicated 
in  Figure  lb.  an  .-Viry  phase  with  very  slow  group  velocity  is  developed,  whicli  is  associa'cd 
with  the  minimum  on  the  group  velocity  curve  (Figure  la),  hor  signals  with  such  a  long 
duic.i ion.  the  spectral  metliod  for  attmuial ion  esiimation  is  not  aiiplicabk'.  1  he  com]CiU' 
v.a'v'e  spectrum  may  not  brr  recovered  by  either  truncating  the  signal  or  taking  the  complete 
v-ave  train.  The  former  ajjproach  removes  a  portion  of  the  wave  energy,  while  the  latter, 
wn.en  ntcasuring  a  sample  of  finite  size,  may  include  reflections  that  bounce  back  and  forth 
between  the  source  and  receiver.  One  may  also  try  to  estimate  attenuation  by  measuring  the 
amplitude  decay  of  the  first  arrival  with  distance.  However,  because  of  the  dispersion  effect, 
tlic  wave  amplitude  decreases  with  distance  even  in  the  absence  of  intrinsic  attenuation  as  in 
this  example  (Figure  lb).  By  studying  this  theoretical  example,  it  is  clear  that  any  attempt 
to  measure  intrinsic  attenuation  using  cylindrical  bars  will  have  to  consider  the  change  of 
waveform  due  to  the  dispersive  nature  of  the  waveguide.  Fortunately,  since  this  effect  is  well 
governed  by  Eq.  (1),  it  can  be  accurately  corrected  for  provided  the  parameters  Vp,  I  ,,  and 
a  of  the  bar  are  given. 

ESTIMATING  ATTENUATION  THROUGH  WAVEFORM 

INVERSION 

In  the  presence  of  intrinsic  at  tenuation  in  the  bar,  the  attenuation  effect  is  taken  into  account 
by  making  the  wavennniber  h  complex,  as 

^•—♦^■(>3-;^)  •  til 
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\v!'..-n'  Q,  is  iho  cxtriisional  (luality  factor  of  th<‘  bar  material.  As  will  be  M  I..-,  i 

ii:  the  .•.\')(nimental  iHoci.Hliire.  the  attenuation  measurement^  am  iisuall_\  uerfu;;!..  !  •. 

ii.'iiKAv  ire()ueni'\'  ratine,  lln'refore.  oven  if  can  vary  with  fn  iiU'r.cy.  it  may  be  :  i  ;e,e.l 
a--  i  iawtaiit  mer  the  narmw  fr('<iu(;iuy  ranj;('.  In  addition,  the  aradastit  wave  disyei • -i 
iM't  included  because  the  clfect  is  negligibh*  in  this  narrow  fierjuency  range. 

I  he  wa\’eform  inrersion  technic|uc  involves  compaiing  the  '.'.areforms  recei\'ed  nn.e  ".'  o 
b.  is  of  tlu'  same  material  but  clifTerenl  lengtlis  .C]  and  .r^  (>  .fj  ).  In  the  two  ■m;.-.  ■;:i- 
n.easurenu'ut.  conditions  ^i.e.,  the  signal  g«‘neiation  and  r^•cei\•e|■  ri'sponse)  are  a^.'ima'i;  O) 
bi  tlie  same.  I'roni  Ibis.  (2)  anti  (3),  it  is  rivulily  seen  that  tin'  ‘spectra  of  the  'wo 
.‘.le  j'elated  via 

lb’(u,',  .rj)  =  ii)  exp(r7:(l  +  ?72C^,)(.r2  -  .r, )]  .  }i 

oi.  transform  into  time  domain 

\V(t,X2)  =  \Vit,Xr)*Dit)*A{t)  .  0  1 

where 

D(t)  =  F~'  {exp[ik{x2  -  Xi))}  and  A{t)  =  F"’ {e.xp[- A-(i2  -  x,)/2Qe]} 

arc  respectively  designated  as  the  dispersion  operator  and  attenuation  operator,  the  symbol 
+  denotes  convolution,  and  •  •}  denotes  taking  the  inverse  Fourier  transform.  From 

Fqs.  (4)  and  (5),  it  is  clearly  seen  that  the  difference  between  the  waveforms  W^t.x^)  and 
ir(<..ri)  is  due  to  two  effects.  The  first  is  the  waveform  distortion  due  to  the  velocity 
dispersion  [/2(t)]  along  ihe.rj-.ri  section  of  the  cylindrical  bar.  The  second  is  the  amplitude 
de<  ay  tine  to  the  intrinsic  damping  [/l(0l  •"'long  the  same  section.  The  first  effect  can  be 
corrected  using  the  theory  presented  in  the  previous  section.  The  second  effect  will  be  ':=ed 
to  derive  of  the  bar  material.  Although  the  waves  in  a  cylindrical  bar  usually  have  long 
durations,  the  attenuation  affects  the  first  few  cycles  of  the  waveforms  in  much  the  same 


way  as  it  affects  llii*  whole  wa\i'  tiain.  Ifascfl  on  lliis  fact,  the  wave  form  inversion  tecliniejui- 
ih  rivcs  attenuation  or  from  tlie  first  few  <  v<  l<  s  of  tlie  wave'form-^. 

I  he  waveform  invcision  pioce<!urc  consists  of  two  maji;r  slejis.  Iji  tJie  lirst.  the;  n,i  a 
sure  waveform  n''(f..ri)  is  tin  orctically  itropagated  to  distance  r-i.  in  which  the  elfcct  of 
.-.t  ti'iiuaiion  is  not  included.  Mathematically,  this  operation  is  cximcssed  as 

ll’(/..iw)  ^  n'(/..j:,)  W;(/)  .  ^(,) 


v...<u<'  11  (f.X;)  is  the  lestiltim.;  waveform  after  the  propagation.  In  this  operation,  if  the 
ities  \  p  and  \  ,  are  o.xact,  \V{t..r2)  will  be  aligned  in  phase  with  the  measured  waveform 
11  (/.  cj;.  In  practice,  a  slight  shift  of  W{t,X2)  may  be  necessary  to  obtain  such  an  alignment, 
because  the  Vp  and  V,  used  in  calculating  Eq.  (6)  are  measured  values  and  may  contain  errors, 
rius  first  step  may  be  called  the  dispersion  correction.  The  second  step  is  the  inversion  for 
attenuation  or  by  minimizing  the  difference  betwc<-n  \V{t,X2)  and  U-'(f,Xj).  To  do  so.  we 
construct  the  following  error  functicjn 


/■To+at  _ 

EiQe)  =  /  [IV'(<,X2)  - 

JTq 

c  To  +  A  r  ,  __ 

=  -  Ii)/2C?,]}]  ,  (7) 


wlierc  To  may  be  chosen  as  the  beginning  time  of  IE(<,X2)  and  AT  is  the  duration  of  tlie 
time  section  in  which  ]V[t,x2)  \V[t,X2)  are  matched.  In  the  inversion  using  synthetic 

Waveforms.  A7’  can  include  any  number  of  cycles  of  the  waveforms  without  altering  the 
inverted  value  of  Q^.  because  the  waveform  data  are  exact.  For  laboratory  data,  however, 
tlie  measured  waveforms  always  contain  experimental  errors  (e.g..  random  electrical  noiseb 
Ihercfore.  AT  should  be  chosen  as  long  as  possilile  to  reduce  the  effect  due  to  the  errors 
In  Jicc.  A7  may  be  chosen  as  2.ri/Fe,  which,  for  the  shorter  samjdc.  is  api)ro.\iinal-  I\ 
trie  arrival  time  difference  between  the  direct  arrival  and  the  first  reflection  back  from  tiie 
source.  Ihc  minimization  is  perfmmed  using  the  non-linear  least  squares  procedure.  We 
fiist  assign  a  very  rough  estimate  of  Then  we  multiply  the  s])ettrum  lV(,.c',  .r ,)  with 
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,  -  .v,)/-XK]  and  tiansfonn  the  product  back  lo  the  time  domain  [because  tid-' 

ration  niainlvM.odiiies.heanuditudr  of  1K(/.X2)  without  (significantly)  altenng  it.  pha  =  e. 

I  li<  -  .  ![<•<(-;  <,r  imdt  iplo  11  llc/l  ions  will  not  be  sliift<'d  into  t  !i<-  intci  \al  A  /  ].  If  the  >  ^t!.l.,•^  cd 
!.),  is  r.ot  sullicieiil  to  make  tl:e  resulting  waveform  m.itdi  with  11  (b.ip).  Q,  :>  ja-i  a, 

fodowim:  the  method  of  non  linear  least  sqnan-s  minimization  (Mo,/..  ]97S)  and  th.-  same 
;  ardinr,'  is  repeat, 'd  in  an  iterative  manner  until  tlic  (uioi  leaches  a  minim.’, m.  1  ..e 

■.,.lue  of  t.b  at  this  minimum  is  taken  as  tlie  estimated  Q.. 

to  summarize  (lie  procedure,  we  give  an  inversion  example  using  synthetic  waveforms. 

1  lie  waveforms  are  generated  using  the  same  parameters  used  in  Figure  1.  But  now  a  0, 

.-30  is  used  in  (lie  calculation  of  the  synthetic  waveforms.  The  waveforms  at  r,  =  lO 
.■m  and  -  -30  cm  are  used  for  the  inversion.  Figure  2a  shows  the  waveforms.  1  he 
waveform  at  10  cm  is  theoretically  propagated  to  30  cm  without  including  the  attenuation 
(i  e.,  Q,  T-  oc).  Figure  2b  shows  that,  after  the  propagation,  waveform  lF(f.T2)  is  aligned  in 
phasi'  with  \V{t..V2).  However,  the  amplitudes  of  the  two  waveforms  do  not  match  because 
(he  attenuation  rlTcct  was  not  included  in  the  propagation.  The  next  step  is  to  minimize  the 
amplitude  difference  between  the  two  waveforms  using  the  inversion  technique.  As  shown  in 
Figure  2b,  only  the  first  two  cycles  are  used  in  the  inversion.  After  the  inversion,  the  two 
waveforms  coincide  with  each  other,  as  shown  in  Figure  2c.  The  Qe  value  obtained  from  the 
inversion  is  50.1,  in  close  agreement  with  the  true  value  of  50.  It  is  noted  that  for  synthetic 
data,  the  inverted  Q,  is  not  sensitive  to  the  number  of  cycles  used  in  the  inversion,  because 
the  data  do  not  contain  errors  or  effects  of  extra  arrivals,  such  as  multiple  leflcctiuns. 

The  pror^yiiirc  of  matching  tlie  first  part  of  n>{/..T2)  and  1V(/,  j^)  using  tile  mvoision 
t.  rhiiiqiie  has  several  major  advantages.  The  first  is  that  only  the  first  few  cycles  r.f  the 
U.iw  furms  .ii,'  nr, 'dr,!  f<u  th,'  Q,  , ’sl imat h-u.  This  sepaiale^  the  effevt  of  inliinsic  an,  n  iamai 
f:om  tin-  effect  due  to  lafer  arrivals  (i.e..  multiple  reiiections).  The  second  advantage  ;h;,f 
the  solution  to  the  inverse  problem  is  uimiue  because  only  one  parameter  Q,  is  estimated. 
Finally  and  most  imiiortantly.  the  iiiveision  is  performed  using  the  waveforms  of  the  low 
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1; ( cmency  fundainciital  UKuie.  allowing  tho  atlcmnilion  to  bo  incasuicd  in  a  nnich  Ic^'aci 
fioquoncy  range  than  that  nf  tlie  ultrasonic  pnKo  ])io])agat ion  technique. 

APPLICATION  TO  LABORATORY  MEASUREMENTS 

1:  th;.'  'Ci.'tion.  we  illn.'-trat .  the  .•tjiplicatii.iu  <>{  the  v.a\'<-furin  iiix'eisiuii  ;ji(;cedine  to  li.e- 

!.;in,iati,ry  nicasureinenl  of  altt.'nualion  using  cylindrical  hnr.s.  The  samples  used  are  a  l’\  (' 
uialeiial  and  Sieura  White  granite.  The  first  is  a  low  velocity  and  highly  attennati\’e  iila-tic 
loalerial.  The  second  is  a  high-Q  rock.  'I'lie  elastic  properties  of  tlie  two  materials  are  given 
in  Table  1. 

Experimental  Procedure 

f  igure  3  shows  a  diagram  of  the  measuring  system.  A  pulse  generator  with  various  source 
functions  is  used  to  apply  an  c.xcitation  signal  to  the  source  transducer.  Unlike  most  pulsi* 
transmission  measurements  where  a  sharp  source  pulse  is  usually  used,  the  present  measure¬ 
ment  uses  a  burst-sine  wave  source  with  adjustable  frequencies.  This  is  because  the  waveform 
inversion  technique  can  deal  with  signals  of  long  duration,  while  other  measurements  (e.g.. 
spectral  ratio  method)  require  signals  of  short  duration.  The  signal  generator  in  Figure  3 
:s  adjusted  to  generate  signals  with  appropriate  frequencies  within  the  frequency  range  of 
the  lundamental  w'ave  mode  in  the  cylindrical  sample.  This  requires  that  the  source  and  re- 
cti\  t’r  transducers  have  good  responses  in  the  (low)  frequency  range.  The  measurements  are 
[ue  formed  on  two  samples  of  the  same  material  but  different  length.  The  requirement  of  the 
•saiire  material  for  the  two  measurements  is  based  on  the  consideration  that  the  transmitting 
a:.>;  n’lxTing  of  the  sourcu'  and  receiver  transducers  depend  on  the  elastic  i)roporties  of  the 
sample,  ''.hich  determine  the  load  impedance  of  the  transdtirers.  Using  samples  of  the  same 
tnateiial  ensures  that  measurement  conditions  such  as  transducer-sampk'  coupling,  sigu.d 
geueiation.  and  signal  receiving  are  tlie  same  for  lioth  measurements,  so  that  the  dilfeieiKo 
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I ■<  t I  II  the  (wo  iix-asuK'd  wa\cforms  nr<-  due  to  tiic  <li^i)cisioii  aiui  attciiuation  effect?  oi..y. 
In  ,ii!ditit»M,  1 1('(  aiisc  (lie  \\a\('torni  inversion  teeliniriue  iec|uue.s  that  llio  two  w;;\<ii,i  n  a- 
li  .Ma  d  in  phase  Iu  IoK'  the  ni\eisiou  is  applii'd.  tlie  sampling  inlei\’al  of  tl;e  digital  0‘^ci.;o- 
:  .  ope  ,-liould  he  fine  enoui'h  to  achieve'  such  an  alignment.  .\  sainjding  rate  o!  aijont  ■ill-lOO 
;  Mil',  pel  c\  lie  is  1(11  MiMi'.elideil. 


!?  ('siilt.s 

!  MW  pi(-ent  the  results  for  the  '^VC  bar.  1  he  lengths  of  the  short  and  long  sample;  are 
'  .‘.r.d  'J1.S2  cm  respixtively.  The  radius  of  the  sample.?  i.s  0.645  cm.  Figure  4a  shous  the 

•  isiir,  d  wa\eforms  for  i  he  S.2S  and  21.S2  em  samples,  d'he  waveforms  are  measured  for  the 
hi  kll,'  and  60  kllz  .source  signal  frequencies.  For  the  shorter  sample,  the  received  (40  kHz) 

■  iL’iial  clearly  shows  the  arrival  of  the  first  and  second  reflections.  The  time  interval  between 
-igual  onset  and  the  arrival  of  the  first  reflection  determines  A7’  in  Eq.  (7).  For  the  longer 
sanqile,  the  received  waveforms  (particularly  the  60  kHz  signal)  exhibit  significant  distortion 
due  to  dispersion  and  amplitude  attenuation  due  to  intrinsic  damping.  For  the  given  material 
inoperties  (Table  1)  and  bar  radius,  the  waveforms  at  8.28  cm  are  theoretically  propagated 
to  21.82  cm  to  correct  for  the  dispersion  effect.  The  waveform  inversion  is  then  applied  to 
derive  from  the  40  kHz  and  60  kHz  waveform  pairs.  The  results  of  waveform  inversion 
are  given  in  Figure  4b.  After  the  dispersion  correction  and  inversion,  the  waveforms  (both 
40  and  60  kHz)  at  21.82  cm  are  satisfactorily  recovered  from  the  waveforms  at  8. 28  cm. 
( 'on.sidcring  the  significant  distortion  of  waveforms  at  21.82  cm  compared  with  those  at  8.2(8 
cm.  tlic  results  in  Figure  4b  show  that  the  dispersion  correction  is  effective  and  suffi-rien* Iv 
accurate.  Furthermore,  the  Qe  \’alucs  obtained  from  the  two  different  frequencies  are 
cl'j'^e.  They  are  Qr  —  14.2  for  40  kllz  aiul  14.9  for  CO  kHz.  indicating  the  consistency  of  tlu- 
iu\<  rsion  rc.sults. 

Of  the  most  interest  are  the  results  from  the  measurements  on  granite.  Resonant  I'ar 

f) 


iii'.'aoUl (  incuts  were  pciforincd  on  the  sample  l>«-)<jre  il  was  i  til  into  two  sainjrles  of  leiigtlis 
5.1G  and  11.91  cm.  I  Ik'  radius  of  llu*  sample.s  is  t).()7‘2  cm.  1  he  resonant  l.nir  results  will  iie 
compared  with  those  from  tlie  ])resent  techni<iue.  Figure  o  sliows  tlie  waveform  data  from 
na  asuriiig  the  shorter  and  longer  samph's  for  dilfcrcnl  source  frccjuencies  ranging  from  id 
to  I-K)  kll;'.  With  incK-'asing  fn'fpieiu  y.  the  dispersive  features  of  the  waxes  become  moie 
(  '.■;do;;i.  l  iguri'  G  shows  the  m.'Gii  of  the  xxavvforms  after  tin'  iiixcrsioii.  The  niati  hes  ,,ie 
gi  ie  i.d!;.'  •.•ery  good.  I  he  iuxcitc'*!  valiK.s  lor  these  frefpieiu  ies  .ii<'  gix'eii  in  Figiiie  o 
for  (..oil  pair  of  niatcheii  waxt  forms.  liar  vahu's  are  on  llu.-  oiihu  (xf  110.  I'o  eoinp.m' 
tie  sc  d.ita  with  those  from  tlie  resonant  har  measurement,  we  plot  the  attenuation  data 
c.\;nessed  as  lOOO/Qf  versus  frerpiency  in  Figure  7,  together  with  a  xalue  measured  on 
a  granite  sample  using  spectral  ratio  method  around  800  kHz.  'I  he  ultrasonic  cxtensional 
attenuation  xalue  xxas  derived  from  tlic  measured  cojnprcssional  and  shear  attenuation  xalues. 
.As  sc'-n  in  this  figure,  the  xvax-eforrn  inversion  results  (open  circles)  are  quite  consistent  xvith 
the  resonant  bar  results  (triangles).  The  ultraisonic  measurement  (square)  yields  significantly 
higher  attenuation  than  the  txvo  low  frequency  measurements,  suggesting  that  scattering 
is  the  mechanism  for  attenuation  at  ultrasonic  frequencies  (Winkler,  1983;  Blair,  1990). 
1  his  may  also  explain  why  the  xvaveform  inversion  results  generally  shoxv  slightly  higher 
attenuation  than  those  of  the  resonant  bar  measurement,  assuming  that  at  the  100  kHz 
fi<  quency  range  tliere  are  still  some  scattering  cfTcrts. 

CONCLUSIONS 

1  hiS  ^t'.idy  jrresents  a  new  apiiroach  to  attenuation  measurements  through  inversion  of  \vax('- 
I'  ';i  ■  In  fact,  the  xvaxeforiii  inxer^ion  technique  is  not  restricted  to  bar  mcasurenionf s.  It 
■  .in  .11  be  .vdapted  to  other  pulse  tiansmission  measurements.  For  example,  for  nltrasc'ni. 
II. ( .iMiK  incnts  u.sing  bam|)le.s  of  large  lateral  ditnensioii.s,  the  waxeform  inxersion  is  apj)li(  .iblc 
if  the  w.;\(  fiirms  are  corrected  lor  difiraction  (or  beam  spreading)  elfects  of  the  transdiuci 
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1  \t  !i  fv  I  tli<‘  1).!!  .i^coiiu't ly.  tlio  tf'<;l)ni<iue  still  has  several  applications.  It  can  be 
..c.ipt,,]  ti.  ii.iasnie  attenuation  unch'r  ('oiifining  pressures  a.j^proiniate  to  the 
■  ‘I  il:o  liopei-ioii  (  liar.u  tei ist i( s  of  tlie  har  in  the  presence  of  a  confinini;  ;nec;;;ni 
1  ■  •  >  I-  I  >'iici  il\  .Kcounletl  for.  In  addition,  the  ua\eforms  of  fleN\iral  ana!  to;oo:;al 
ei  .(  I  ■,  in.iii  it  al  bar  can  also  be  utilized  to  estimate  the  shear  attenuation  of  ti:c 
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laM('  1:  Density  /),  eoinpressional  velocity  \  p,  and  shear  velocity  of  the  solid  ii-cd  in  the 
iiu'asiireuu'ut. 
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FREQUENCY  (kHz) 


b. 


('  1;  (aj  Phase  an<i  Rioiip  velocities  ((iaslied  ciuves)  "f  tlio  cyliucliical  har  (u  1  i  ;a. 
1  P  ---IJ72  iii/s,  aiui  \',  -'J. ')()()  Ill/s)  in  tie'  frefm<'n<:y  ra !!{;<■  t:o\'ei ed  l^y  t li<‘  soui i.e  s]iect udd 
(solid  curve).  I  he  vclociU<s  are  sealed  l)y  m/f.  (b)  Proi)agation  of  the 

fundamental  mode  in  the  bar.  Note  the  distortion  of  waveforms  and  the  dev'clopment 
of  the  Airy  phase. 


a . 


so  70  90  '.10  130  ISO 

TIME  (WicroSec) 


I'  ij'.iiK’  2;  lllustialioii  nf  the  wiivt.'foiiii  inversion  procedure.  (a)Piopagate  waveform  at  jj  to 
.cj.  (l?)Rcsult  of  the  ilisiicrsion  correction;  the  waveforms  are  aligned  in  phase;  the 
amplitude  difTciencc  is  due  to  attenuation  and  is  to  be  minimized  to  find  Qt-  (c)  After 
inversion,  the  two  waveforms  are  matched.  The  required  for  the  match  gives  the 
estimated  attenuation. 
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SAMl’LK  (cvlindricijl  bar) 


SOURCK  TRANSDL’CKR  RECEIVER  TRANSDUCER 


F;gni("  'U  Diagram  of  tho  system  for  measuring  attenuation  in  bars. 
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n. 


Measured  waveforms  from  PVC  bar  (a=0. 645cm) 


b.  Results  of  waveform  inversion 


igurc  1;  (a)  \Vav(.'forms  iiu-asiirccl  from  (he  short  and  long  P\'C  bars  at  -10  and  CO  kHz 
frequencies.  (b)Ma(thcd  waveforms  after  tlie  dispersion  correction  and  inversion.  The 
iinerlcil  values  for  the  two  freriurncics  are  also  indicatefl. 
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Waveform  Inversion  Results  for  Sierra  White  Granite 


0  20  40  60  80 

TIME  (MicroSec) 


1110  6:  Matched  wavchaiiis  after  tlio  <lis])orsion  correction  and  inversion  for  the  diffeient 
frcciuencics.  I  he  inverted  Qe  values  for  the  frequencies  are  also  indicated. 
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ATTENUATION  (1000/Q) 


COMPARISON  OF  ATTENUATION  VALUES  OF  SIERRA 
WHITE  GRANITE  OBTAINED  USING  DIFFERENT  METHODS 


5  0  r 

4  0  I- 
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ultrasonic  pulse 


v,raveform  inversion 


10  100  1000 

FREQUENCY  (kHz) 


r  i:  Couip.'ii  ison  of  extol,  ■i’jnal  attemialiuns  (exiiicsse'l  as  1000/Q)  irica?ui(.-d 
I'l-uiianl  bar  (tnanglos J.  \\a\t  frain  inversion  (open  rircks),  and  ultrasonic  pulse  jirop- 
a  gat  ion  (square). 
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